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Abstract 

We consider the totally asymmetric simple exclusion process, a 
model in the KPZ universality class. We focus on the fluctuations of 
particle positions starting with certain deterministic initial conditions. 
For large time t, one has regions with constant and linearly decreasing 
density. The fluctuations on these two regions are given by the Airyi 
and Airy2 processes, whose one-point distributions are the GOE and 
GUE Tracy- Widom distributions of random matrix theory. In this 
paper we analyze the transition region between these two regimes 
and obtain the transition process. Its one-point distribution is a new 
interpolation between GOE and GUE edge distributions. 



1 Introduction 

In the search of universal hmit distribution functions and hmit processes, 
we consider the KPZ universahty class (KPZ for Kardar-Parisi-Zhang) orig- 
inally introduced for stochastic growth models [13]. For growth in 1 + 1 
dimensions the scaling exponents of fluctuations, 1/3, and correlations, 2/3, 
can be (non rigorously) determined by some involved arguments, see e.g. [16] 
for an extended discussion. However, to get more insights into the limit laws 
and limit processes, one is led to consider solvable models in the universality 
class. 
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Figure 1: Left: the density p for large time t is linearly decreasing from 
(0,1/2) to (t,0). Right: The limit shape in an associated growth, obtained 
from the density (/c G N is the label of the particle which starts at —2k and 
Xfc(t) is its position at time t). 



One such model is the polynuclear growth (PNG) model in discrete time, 
which has two interesting limits, where new processes have been discovered. 
The first limit is the continuous time PNG model, for which it has been 
shown that the surface growing in a droplet shape is, in the large time limit, 
governed by the Airy2 process [20] . The second one is the totally asymmetric 
simple exclusion process (TASEP), in which quite recently the limit process 
of the particles positions starting from a periodic initial conditions has been 
unravelled and called Airyi process [2,21]. In the surface growth picture this 
corresponds to the flat initial conditions. 

If Xk(t) denotes the position of particle with label k, one of the usual 
geometric representation of the TASEP in terms of surface growth is obtained 
by the graph {{k,Xk(t) + k)}, see Figure [T] (right) and also e.g. [5]. By 
universality it is expected that the limit process in one-dimensional KPZ 
growth is the Airy2 process for the curved regions of the limit shape, and the 
Airyi process for the flat parts. However initial conditions can easily generate 
limit shapes which have both curved and fiat regions. Therefore there exist 
transition regions where the limit shape smoothly changes between curved 
and fiat. 

The novelty of this paper is the analysis of this transition region in the 
framework of the TASEP. The observables we consider are positions of several 
particles at time t. In [10], step- initial conditions (particles starting from 
Z_) have been considered from the perspective of a growth model and it was 
proved that the Airy2 process appears in the large time limit. In [1,2] we 
considered periodic initial conditions (particles starting from ciZ, d = 2,3, . . .) 
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and obtained the Airyi process as the hmit process. 

To obtain both regimes and the transition region, we consider in this 
paper particles starting from 2Z_ as in [21]. There are four regions of interest 
as illustrated in Figure [H (left). 

(1) Constant density region. The limit process of particle positions is 
given by the Airyi process Ai- In particular, the one-point distribution is 
Fi(2^/^s), with Fi being the GOE Tracy- Widom distribution. 

(2) Linearly decreasing density region. The limit process is the Airy2 
process A2, which has F2{s) as one-point distribution, with F2 being the 
GUE Tracy- Widom distribution. 

(3) Finite distance from the right-most particle. There the particle 
positions are described via the GUE-minor kernel [12]. In particular, the n-th 
right-most particle is distributed as the largest eigenvalue of the n-particle 
GUE ensemble. 

(4) The transition region between (1) and (2). The fluctuations are 
governed by a new process obtained in this paper: the transition process 
Airy2^i, which we denote ^2^1- In particular, the one-point distribution 
interpolates between F2{s) and Fi(2^/^s) and the transition region has width 
which scales in time as t'^^^. 

The analysis is done by using the framework of signed determinantal point 
processes introduced in [2]. This new approach allows us to analyze all four 
regions for our initial conditions. This is contrasted to the previously used 
determinantal point process issued by the RSK construction, by which only 
the step initial condition or its variants could be analyzed [3,4,6,8-10,20,22]. 
We explain how the analysis has to be done for all 4 cases, but the complete 
asymptotic analysis is presented only for the transition region, the technically 
most difficult one, and the really new result of this paper. The result is a 
process, A2^i, interpolating between the Airy2 and the Airyi processes. For 
more details about the Airy processes, see the review [5]. 

The transition we discovered is not the first one between some GUE 
and GOE type distributions, but it seems to be different from the one pre- 
viously known for random matrices, non-colliding Brownian motions with 
open boundary condition and so on [7,14,19,22]. The main differences are 
the following. On the natural scale of the problems considered, the final dis- 
tribution is Fi(2^/^s) for our case and Fi{s) in the previous case. Secondly, in 
the previous case, the GOE- type distribution appears at a single point, while 
in our case, the GOE-type distribution is on an extended region. Moreover, 
our transition smoothly interpolates between F2{s) and Fi(2^/^s), which is 
not the case for the other transition. In principle, we can not however yet 
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exclude that by a change of variable, with both the rescaling of fluctuations 
and spatial correlations, the two transitions map one to the other. 

Outline. In Section[2]we define the model we analyze and state the results. 
In Section [3] we explain the finite time result and set the scaling limit. In 
Section H] we do the complete asymptotic analysis for the transition region 
and in Section^lwe explain how to do the analysis for the other cases. Finally, 
we present an explicit form of the transition kernel in terms of Airy functions 
in Appendix |A] and we explain the correctness of the Fredholm determinants 
involved inlBl 
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2 Model and results 

In this paper we consider the continuous-time totally asymmetric simple 
exclusion process (TASEP) on At any given time t, every site j & 1i can 
be occupied at most by one particle. Thus a configuration of the TASEP can 
be described hj r] = {rij,j G Z|r7j G {0, 1}} G = {0, 1}^. rjj is called the 
occupation variable of site j, which is defined by rjj = 1 if site j is occupied 
and rjj = if site j is empty. 

The dynamics of the TASEP is defined as follows. Particles jump on 
the neighboring right site with rate 1 provided that the site is empty. This 
means that jumps are independent of each other and are performed after an 
exponential waiting time with mean 1. More precisely, let /: f2 — R be 
a function depending only on a finite number of r^^'s. Then the backward 
generator of the TASEP is given by 

Lfiv) = E^^(l - VHi){fiv"^') - m)- (2.1) 

Here 7]^'^~^^ denotes the configuration rj with the occupations at sites j and 
j + 1 interchanged. The semigroup e^* is well-defined as acting on bounded 
and continuous functions on Q. e^* is the transition probability of the 
TASEP [18]. 
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We denote by Xk{t) the position of the particle number k at time t. As 
initial condition we consider particles starting from 2Z_, i.e., Xk{0) = —2k 
for k = 1, 2, . . .. On the macroscopic level, the limit particle density ■u(^) is 
given by 

d 1 f 1/2, e<0, 

«(0 = ;T7lim-E(#(fc:a;,(t)>et)) = <^ 1/2-^/2, ^€[0,1], (2.2) 

dO-oot 

Thus for large time t the expected number of particle at site x, p{x), is close 
to u{x/t), see Figure [TJ 

As observables we consider the positions of finite subsets of particles, 
{xi{t),i G /} for some / C N, |/| < oo. The scaling limits we have to take 
depend on which of the four regions described in the Introduction we focus 
on, see also Figure [H The main result of this paper is the description of the 
large time fluctuations in the transition region, which now we describe. 



(4) Transition region: the A2^i process 

The transition region has width of order t^/'^, which is indicated by the fact 
that the index of the particles which at time t are around x = fluctuates 
on the t^/^ scale around the macroscopic value t/A. Therefore we set 

niT,t) = [t/A + Tit/2r/% (2.3) 



The density f l2.2p changes in the transition region. The limit density can be 
used to determine, on the macroscopic scale, the expected location at time t 
of a particle with index n{a) = [t/4 + at] (of course a > —1/4). The result is 
then 

lim ^ = / ^ - + « ^ [-1/4, 0], ^2 4) 

t^oo t \ -2a, a > 0. ^ ' ^ 

By using (12.41) with at = r(t/2)^/^ we are led to define the rescaled process 
of particle positions by 

^ , , _ Xnir,t){t) - (-2r(t/2)^/3 + min{0,r}^(t/2)V^) 
^ ^ Mr) - _(t/2)i/3 • ^^-^^ 

The main result of this paper is the convergence of Xti^r) to the transition 
process ^2^1 defined below. 

Definition 1 (The Airy2^i process). Let us set 

s- r- > 

s.-ri, n<0. (2.6) 
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Figure 2: An illustration of the paths 7-|_ and 7_ in Definition [H 



and define the transition kernel 

KociTi,si;T2,S2) = , , , exp -— ^ ^ 1{t2 > Ti) 

^/4n{T2 - Ti) V 4(r2-ri)y 

wzt/i t/ie pat/is 7+,7- satisfying —7+ C 7_ wi/i 7+ : e^'^+oo — >• e~"'^+oo, 
7_ : e~"^"cx3 — * e^'^^oo for some 0+ G (7r/3,7r/2), 0_ G (7r/2,7r — 0+), see 
Figure 

T/ie Airy2_+i process, ^2-*i; process with m-point joint distributions 
at Ti < T2 < . . . < Tm given by the Fredholm determinant 



Pf P|{^2^l(rfc) < Sk}j = det(l - XsKooXs)LH{Tl,...,Trr 



}xR) 



(2i 



where Xs{Tk,x) = l{x > Sk). An explicit expression for K^o in terms of Airy 
functions can be found in AppendixlAl 

Remarks: A2^iit + r) becomes 2~^/^^i(2^/^r) as t — cx) and ^2(t) when 
t —00. The Fredholm determinant in (12.81) is well defined because, as 
proven in Proposition [9] of Appendix [B], there exists a conjugate kernel of 
XsKooXs which is trace-class on H = L"^ {{ti, Tm} x R). 
Now we can state precisely our main Theorem. 

Theorem 2. The convergence of Xt to the transition process A2^i, 

limXi(r) =^2^i(r), (2.9) 

t~*oo 

holds in the sense of finite- dimensional distributions. 
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A remark on initial conditions. In this work as well as in many of the 
previous papers in the field, the situations analyzed with deterministic ini- 
tial conditions might look quite pecuhar: step- initial conditions [11], periodic 
with period 2 or more [1,2]. However, it is intuitively clear that small pertur- 
bations of the initial conditions do not affect the large time behavior. This 
is indeed the case by a coupling argument. 

Consider two TASEP initial conditions of N particles, A(0) = {xj\f{0) < 
... < X2(0) < a;i(0)}, Z{0) = {zn{0) < ... < ^2(0) < ^i(O)} with 
A(0) < Z{0) meaning Xk{0) < Zk{0), k = 1, . . . , N. By a standard cou- 
pling argument, see e.g. [17], for any subset / C {1, . . . , N}, 

P{{xi{t) <a,,ie I}) > P{{z^{t) <a„ie /}). (2.10) 

We can apply fl2.10p to our case to show that the limit result is unchanged 
if we do any bounded perturbation of the initial condition. In Theorem [2] 
we started with initial conditions Xi{0) = —2i. Consider any other initial 
condition Z = {^^(O)} and define 

M = max{\xi{0) - Zi{0)\}. (2.11) 

Then, by fl230D . we have 

P{{x,{t) <ai + M,te I}) > P{{zi{t) <a,,te I}) 

> Piixiit) <ai-M,iEl}). (2.12) 

In the scaling limit fl2.5p . the first and last term in fl2.12p have the same limit 
as t ^ 00 as long as limf^oo M/t^/^ = 0. This holds in particular if Z is any 
bounded perturbation of X, i.e., if M < cx) is independent of t. 

For completeness we state the results in the other three regions. In Sec- 
tion [5] we outline how the asymptotic analysis for the transition region has 
to be modified in order to obtain the results. The scaling is obtained using 

(1) Fixed particle number: GUE(n) minors 
Consider particles with index not rescaled in time, i.e., 

n of order one, (2.13) 
and the rescaled random variables 

Xt{n) = i^^^. (2.14) 



Then, in the t —>■ 00 limit, one gets the GUE-minors(r;,) given in [12], 



hm Xt{n) = GUE-minors(n). (2.15) 

t— »oo 
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(2) Linearly decreasing density region: Airy2 process, A2 

For < a < 1, define 

n{T,t) = [at/4 + T{t/2f/% (2.16) 

and the rescaled process 

" ^ = ^(^T^F^ • 

Then in the t ^ oo hmit, one gets 

\2/3 



(2.17) 
(2.18) 



hm X,(r) = ^^'' ^,(ra^/^(2 - v^)^^). 

(3) Constant density region: Airyi process, Ai 

For a > 1, 

n{T, t) = [at /A + r{t/2f/^] , (2. 19) 

and the rescaled process variables 

. „ ^, ^ .„,..,W-((l_-a^t/2-2.(V2)V3) ^ ^^^^^^ 

Then in the t oo limit, one gets 

lim Xtir) = 21/3^1(2'/^). (2.21) 

t— >oo 

3 Kernel and its scaling limit 

In this section we derive the expression of the joint distributions of particle 
positions and then set the proper scaling limit. 

Consider particles starting at time t = at positions Xk{0) = —2k, 
k = 1, . . . , N. In Theorem 2.1 of [2] we proved that the joint distribution 
of the positions of the particles are given by a Fredholm determinant. The 
kernel is determined via a certain orthogonalization, which for our initial 
conditions has been made in Lemma 4.1 of [2] (with z = x + 2n~2N replaced 
hj z = X + 2n). Once the orthogonalization is made, one can compute the 
kernel which is (4.11) of [2] (with Zi = Xi + 2ni — 2N replaced by Zi = Xi + 2ni). 
This is summarized in Proposition [31 
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Proposition 3. Let particle with label i start at Xi{0) = —2i, i = 1, . . . , N . 
At time t, the particles are at positions Xj. Let cr(l) < o"(2) < . . . < o"(m) 
he the indices of m out of the N particles. The joint distribution of their 
positions Xcr(k){t) is given by 

m 

P( n {Xa{k){t) > flfc}) = det(l - XaKtXa)e^{^il),...,^im)}xZ) (3.1) 

fe=l 

where Xa{cr{k),x) = l{x < a^). The kernel Kt is given by 

Kt{ni,Xi]n2,X2) = ^ ) l[«2>ni] + ^tl^-i' ^^i; ^2, Xa), (3.2) 

\n2 -ni-lj ' 



where 



Kt{ni,xi]n2,X2) = f dv (j) du — (3.3) 



(2vri) 



ro 



g^t^ni 1 + 2t; 

(1 + it)^i+"i+i (^u-v){l + u + v) 

where Tq, resp. T_i, is any simple loop, anticlockwise oriented, which includes 
the pole at V = 0, resp. u = —1, satisfying — 1 — Fq C r_i, i.e., all the points 
of —1 — Tq lie inside the loop T_i. 

In order to prove Theorem [2], we need to focus at particles with number 
Hi close to t/4 since these particles will be in the transition region at time 
t. The transition region has width which scales as t'^^^. The limit density is 
constant to the left of the transition region and it is decreasing linearly to 
the right of it. Therefore, the scaling limit used to prove the main theorem 
is 

m = [t/4 + n{t/2n 

X, = [-2n{t/2f'-~s,{t/2y/% (3.4) 



(3.5) 



where 

Si, n > 0, 

Si - rf, Ti < 0. 
As a consequence the rescaled kernel writes 

Kr'{ri,s,;T2,S2) = Kt{m,x,;n2,X2)it/2y/^2-'-^^ (3.6) 

where 2^^"^^ is just a conjugation so that the kernel has a proper limit. 
We denote by Kl'^^'^ the term of the rescaled kernel without the binomial 
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contribution, which then writes 



exp(t/o(n) + (t/2)2/3r,/i(n) + (t/2)V35 + /3(^))' ^ " ) 

where the functions /, are given by 

fo{v) = -v + lH{l + v)/v), 

h{v) = -H-Av(l + v)), 

h{v) = -ln(2(l + t;)), 

h{v) = ln{l + v). (3.8) 

From now on the r^'s are some fixed values. With this preparation we can 
proceed to the asymptotic analysis needed to prove Theorem [21 



4 Asymptotic analysis 

Proof of Theoreml^ The proof of Theorem [2] is identical to the one of Theo- 
rem 2.5 in [1], provided the following Propositions HI El El H and El (conver- 
gence on bounded sets and large deviations bounds) hold. □ 

Proposition 4 (Uniform convergence on bounded sets). Fix any L > and 

Xi,Si with the scaling ^3.4^ . Then, uniformly for (si,S2) G [— L,L]^, 



lim /sTf '^(m, xi; na, xa) = K^'^'in, si, T2, S2) (4.1) 

t— >oo 



where 



(4.2) 

with the paths 7+, 7- satisfying —7+ C 7_ with 7+ : e^'^+oo — * e "'^+00, 
7_ : e~^'^-oo e^^^oo for some 0+ G (tt/S, vr/2), 0_ G (7r/2, vr — 0+). 

Proof. The first step is to control the contribution away from the critical 
point given by 

^ = _iW=0 ^ . = -1/2. (4.3) 

dw 4:v{l + v) ' ^ ' 
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If we write v = x + iy, x,y E H, then we can analyze 
Mfoiv) - /o(-l/2)) = -{x + 1/2) + I ln(((l + xY + y')/{x' + y')). (4.4) 

This expression equals zero for 

a) X = -1/2, y eR, 

b) y = ±g{x), with g{x) = -^ii^^^^^f^^. 

If is easy to see that the solutions ±5'(x) are symmetric with respect to 
V = —1/2 and they go around —1 and once. Moreover, the loops leave the 
critical point f = — 1/2 in the directions e"^™^^ and e^^^'^/^, see Figure El We 




denote by Di, . . . , the following regions: Di is the region enclosed by ±g 
around —1, D2 is the rest with real part less than —1/2, D4 is the symmetric 
image w.r.t. —1/2 of Di and D3 of D2, see FigureO Then Fq can be chosen 
to be any simple anticlockwise oriented finite length path staying in D3 and, 
similarly, F_i is chosen to stay in D2 (except at f = —1/2). The constraint 
— 1 — Fo C F_i is easily satisfied except that for Fq we have to go through D4 
too, very close to f = — 1/2. Moreover, we can take Fq leaving from —1/2 
with an angle between —ir/Q and — 7r/3. Similarly, F_i leaves in the direction 
from 27r/3 and Svr/G. This will simplify the argument for moderate and large 
deviations. 

Let us set Tf^ = {v e Fq, |t;+l/2| < 6} and Fl^ = {u E F_i, |m+1/2| < S}. 
Then the integral is over Tq U F_i = Fq U F^^^ + S, where S is the rest of the 
contours. The first step is to bound the integral over E. For < 5 <C 1, we 
can choose Fq and F_i such that, for (m, v) G S, \u — v\/6 and |H-u + f |/5 are 
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Figure 4: The paths Fq and F_i close to the critical point —1/2. The dashed 
lines are the zeros of Re(/o(x + iy) — /o(— 1/2)). 



bounded away from 0. Then, on E, 



l+2v 



{u—v)(l+u+v) 

0{1/S) and, for some cq = co{6) > 0, Re(/o(f — /o(— 1/2)) < — cq and/or 
— Re(/o('u) — /o(— 1/2)) < — Cq. Thus, the integral over S can be bounded as 



<\u — v\ -^ + |l + 'U + f| ^ 



(4.5) 



for some Ci > 0. For t large enough, both e'^^*^'^^) and Cit^^^ are bounded by 
g-cot/4_ 'jj-^ug^ for t large enough, we have the bound 



(4.6) 



The second step is to control the integral over Fq U F1^. Since S is small, 
we can apply Taylor series expansion on the functions fi defined in (13.81) . For 
this we change variables by setting 



M = -l/2 + f/, v = -l/2 + V 

and we denote 7^ = Fg + 1/2, 7! = F^ ^ + 1/2. We have 

/o = i + if + |y^ + 0(\/^), 

/i = W^ + 0{V'), 

h = -2V + 0{V^), 

h = -H2) + 0{V). 



(4.7) 



(4.^ 
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Therefore the integral over Fq U Fl^^ is given by 



X e 



{t/2f/^ f f g|iV3+(t/2)2/3r24F2-S2(t/2)l/32y 

(27ri)2 j^s ^^(f/ ~V)iU + V) e|tt/^ + (V2)=^/«n4l/2-si{t/2)l/32[/ 



(t/2)l/3 f f W g|<V3+(t/2)2/3r24y2-S2(t/2)i/32y 

(27ri)2 ^^(f/ - + V^) e|tt/3 + (V2)2/3ri4C/2_5^(t/2)l/32[/ 



+R. 



(4.9) 



To bound the remainder, R, we use |e^ — 1| < |x|e'^' apphed to x = 
Moreover, note that 0{t'^^^V^) is dominated by 0{tV'^). Therefore, 



1^1 < C2t^/^ / dV dU 



li Jit 



AV 



63 



|ty3+(t/2)2/3r24V2-S2(t/2)i/32y 



_V){U + V) e|tt/3 + (t/2)2/3ri4C/2-si(V2)l/32[/ 



X e 



(4.10) 



At this point we do the change of variables V = w{4:t) and U = z{4:t) 
and obtain 



\R\ < cst"^/^ / dw I dz 

'(4t)i/3^« J(4t)i/3^« 



W 



/3+T2W^ — S2W 



{z - w){z + w) e^^/3+^i^^-*'i^ 



X e 



c(w;^ z^ L;2^ z) 



(4.11) 



By choosing S small enough, we may assume that 0{wH^^^^) <^ w^, 
0{zt~^/^) < 1, and for t large enough 0{Lt~^/^) < 1. Therefore, 
the exponential in the integral in the w variable can be bounded by 
I exp(xow^^/3 + T2XiW7^ — 52X2"^)! for some Xo^Xi^X2- By choosing 5 small 
enough, the x's can be made as close to 1 as desired. More importantly, for 
5 small, one has Xo > 0. Similar for the variable z for some Xk- We have 



1^1 < csr^/^ / dw I dz 



X 0{w\Lw\z\Lz\z) 



w 



{z -w){z + w) -^1X2^ 

(4.12) 



The integral in (14.121] . without the prefactor t is uniformly bounded 
in t. In fact, the only dependence on t is at the boundaries of the integrals, 
which are at 5e^'^+ and 5e^^^~ with 6+ G (7r/6,7r/3) and 6. G (27r/3, 57r/6). 
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The convergence is ensured by the fact that Re(w'^) = cos(36'+), with 
cos(3^+) < 0, and Re{-z^) = -5Hcos{39_), with cos(3^_) > 0. Thus, the 
and 2;^ terms dominate the others at the boundary of the integrals and 
this domination becomes stronger while t increases. The final result is that, 
we can set 5 > small enough and then for t large enough we have 

\R\ < c4~^'^. (4.13) 

The last step is to analyze the first term in r.h.s. of fl4.9p . One does the 
same change of variable as above and gets 

/ ^ dz ,^ , „.3/.^.,.2_,,. • (4.14) 



(27ri)2 7(4^)1/3^^ JiAtyrs^t {z -w){z + w) e-Vs+n.^.,,, 

We can extend the paths to t = 00 and by doing so we gain the error term of 
order (9(e~^^^^*) for some C5 > 0. With this extension the paths satisfy the 
conditions of 7+ and 7_ of the Proposition. 

Just to summarize, the error term we have accumulated during the above 
procedure is 

0{^5-i^-cot/2^ C4ri/3, e-^^^'*). (4.15) 

□ 

Proposition 5 (Moderate deviations). For any L large enough, 3eo(-^) > 
and to{L) > such that, \/ < e < Eq and t >to, the estimate 

Kr{rusv,r2,S2)\ < e-(^^+^^) (4.16) 

holds for (si, S2) E et2/3]2 y ^2 ^ 

Proof. In this proof we introduce the notation, (jj = Sit~'^^^2~^^^ E (0,^], 
i = 1,2. We divide the analysis in the cases §1 > S2 and si < §2- The 
strategy is the following. First, for the case si > §2, we choose the same 
paths Tq and r_i as in Proposition H] except for a small modification close 
tof = u = — 1/2. We then see that in the unmodified part of the paths one 
has the same integral as for the case o"i = o"2 = times a factor which can be 
simply bounded and gives the needed decay. Then we consider the modified 
parts of the integration paths and see that the integral over these has also the 
required decay. Secondly, for the case Si < S2, we first modify the condition 
of the integral since, otherwise, the optimal paths for the exponential can 
not be followed close to the critical points. The modification produces an 
extra ter m, a residue, which is a simple integral and it can be bounded in a 
similar way. 
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Case (Ti < a-i- The paths Fq and r_i as represented in Figure [5l 

The modification with respect to the ones in Proposition H] is just one 

vertical piece, given by Fyert = + v^(l + ^ ^ [— «]} for some 

ae (1/V3,V3). 

With respect to the case a\ = = 0, the integrand in the integral 
representation of the kernel Kf^'^, see (13 .Zp . has the extra factor 

exp(-ta2 ln(2 + 2v)) exp(tai ln(2 + 2m)), (4.17) 

whose magnitude is given by 

|(j3j7])| = exp(-tor2ln(2|l +f|))exp(tcTiln(2|l + M|)). (4.18) 

a) For the term (1 + 2v) / {{u — v){l+u + v)), we can choose t ^ 1 such that 



dist(Fo,-l -F_i) = 




(4.19) 



which is much better than in Proposition 4, where we had, see Figure H] 
dist(Fo,-l -F_i) = 5t~^''^. Therefore the term {l + 2v) / {{u-v){l + u + v)) 
does not create any problems. 

b) Similarly, F_i can be chosen such that the maximum of |l+u|, for u G F_i, 
is obtained at u = —1/2, thus 

gi<nln{2|l+«|) < ^_ (4.20) 
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c) To can be chosen such that the minimum of |1 + f | for i; G Fq \ T^ert, is 
obtained at r^ert for ^ = ±a. A simple computation leads to 

for e small enough. 

d) Now we evaluate the integral over T^^rt- As considered in the ^ variable, 
the prefactor t~^/^ cancels out and ^ varies over an interval of order one. 
Therefore, to estimate the integral it is enough to estimate the integrand. 
Since e is small, a2 is small too. Thus, Fyert is very close to —1/2 and we can 
apply Taylor expansion of the integrand. The term with the exponential in 
the V variable becomes (f = — 1/2 + ^/o2{l + iO/2) 

exp (t/o(-l/2) + |taf (l+iO'+r2(t/2)2/V2(l+iO'-^2^'t(l+iO + (^(i^2))- 

(4.22) 

By using (T2 = S2t^'^^^2~^^^ and computing the real values of the exponent, 
we get 

irai < exp (t/o(-l/2) + 5f 2"V2(_| _ 1^2 ^ ^ _ ^2^y 

(4.23) 

Here we have S2 > L, thus S2 > 52/2 for large L and 3> S2. Therefore, 
the integrand to be studied can be bounded by 



,t/o(-l/2)^-Sp/2 



(4.24) 



for L large enough and e small enough. The factor e*-^"'-"^/^-' is cancelled 
exactly with the one coming from the integrand in the u variable. 

For the case di = (72 = 0, the analysis of Proposition H] leads to the bound 



on the kernel K^"^^^ 



< C6 (4.25) 



for some constant Cg > 0, as soon as t is large enough. 

Putting together the results of a)-d), the kernel is bounded by cq times 
the factor e~^'^ For L large, §2 > L/2 and S2 > S2/"\/2, therefore 

cee"^'''"/' < cge-i^'^ < C6e~i^('^+'^) < e-^''+'^'^ (4.26) 
where we used S2 > si. 
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Case o"i > o"2- To obtain the bound for this case, we use a different 
expression for the kernel KY"^^, namely 

= itl2f'^-\- I dv I du- \^ (4.27) 

* (27ri)2j;„ {^u-v){l + u + v) ^ ' 

exp{tfo{u) + {t/2Y/HMu) + {t/2y/n,h{u) + h{u)) ^ 

where 

^ (t/2)i/3Zi / ^y(,t{h{v)~M~i~v))^(t/2)y2{T2h(v)-Tih{-i~v)) 

with the constraint r_i C — Fq instead of — Fq C F_i. The term I2 comes 
from the fact that, for any fixed f , the new constraint on the paths is obtained 
by deforming F_i and during this process one passes via a simple pole at 
« = — 1 — f , whose residue is I2. 

The analysis of the double integral term in (14.271) is the same as in the 
previous case, where however (m, si, ri) play the role of (f , S2, T2), so this time 
it is F_i which is modified instead of Fq (symmetrically w.r.t. —1/2). We can 
then get as in (14.261) the bound exp(— (si + S2))/2 and it remains to prove 
that I2 is bounded by exp(— (si + S2))/2 too. 

Denote ho{v) = fo{v) - /o(-l - v). It is given by ho{v) = -1 + 2fo{v). 
Therefore the regions where sign of Re{ho{v) — ho{—l/2)) is positive and 
negative are again the ones of Figure [3l In the case di = cr2 = 0, one can 
do essentially the asymptotic analysis made to obtain the estimate on the 
integral over F_i of Proposition H] and we get that the integral is bounded 
in the t — 00 limit. The corrections to the limit expression are of just order 
0{t~^^^ , e"^^), for some /i > 0. But here we are in the case Si G [L,et'^^^]. 
The difference with respect to the case ai = o"2 = is a factor of magnitude 

exp(t(Tiln(2|i;|) -ta2ln(2|l + w|)), (4.29) 

in the integrand. The Fq used for the Ui = cr2 = asymptotic analysis can 
be chosen such that, while going away from the critical point v = —1/2, 

a) \v\ decreases, thus ln(2|f |) decreases, 

b) |1 + v\ increases, thus — ln(2|l + v\) decreases, 
take for example — 1 — F_i of Figured 

Now we use the same trick as above, namely we modify Fq only in the 
neighborhood of 1; = —1/2 as in Figure [5] (just this time the distance to 
V = —1/2 is Y^/2 instead of y^/2). We denote F^ert the vertical piece 
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here too. Then, the contribution on Fq \ T^ert carries an extra term (as in 

g-tailn(2|l+.|) <g-si/^/2^ (4.30) 

for e small enough. Then, for L large enough, | (14.301) | < — e^'^^^(^i+''2) ^^t^ 
some C7 > 0. 

For the contribution of the integral over r^gi-t , we set v — — \ / 1 -\- V and 
do Taylor expansion. Then set V = + i^)/2 with C, G [—a, a], for some 

a G (1/v^, v^). The integral over Fvert is an integral over [—a, a], which 
writes 

X e-K+'^2)v^T(i+iOt{i+ci(v/£))gCi(V^)_ (^4^31) 
We then use 

a) Re((l + iO=^) = l-3e^ 

b) Re((l + iO^) = l-e^ 

c) > ^/L, 

d) 2/3 < |1 + C(v^)| < 2, for e small enough, 
to obtain that | (14.311) | is bounded by 



d^csv^exp (taf (i -e- C9(l - e')/v^)) exp ( - f (^i + a,)t^,). 

(4.32) 

The integral (I4.32p is bounded and, for L large enough, the integrand is 
maximal at ^ = 0. Thus 



<SM) < Ciov^sTexp -C9//L- |((Ti + (T2)v^t) 

< exp(-i((Ti + a2)v/^t) (4.33) 



for L large enough. Reinserting the expressions for cxi and a2, we have 

|(I132DI < exp ( - cio(5i + S2)t^/h) < exp ( - Cn(si + S2)tVL) (4.34) 

for L large enough and some Cn > 0. This bound is good enough to get 
exp(— (si + S2))/2 as bound for L large enough, e small enough and t large 
enough. □ 

Proposition 6 (Large deviations). Set e > 0, then for t large enough we 
have 

Kr"in,s,;T2,S2) <e-(^^+^^) (4.35) 



for (si,S2)G[-L,oo)2\[-L,£t2/3] 



2 
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Proof. One can do large deviations directly, but a shorter way is to use the 
result of the moderate deviations. As in the proof of Proposition we use 
the notation, cij = Sjt~^/^2~^/^, i = 1,2. 

Case (Ti < (72. The term linear in t in the exponential is 
exp(t/o,<x2(^) -tfo,ai{u)), where fo,a{v) = fo{v) - crln(2 + 2v). To obtain 
the bound we just remark that 

fo,M = fo,e/2{v) - {a2 - e/2) ln(2 + 2v). (4.36) 

We take Tq to be the one used for moderate deviations with cr2 = e/2. Tq 
satisfies |1 + f | > 1/2+ a/^/2/2. o"2 > £ implies a2 — e/2 > (J2/2. Therefore, 
for e small enough and t large enough, 

\exp{-t{a2-e/2) ln(2 + 2t;))| < exp(-it(T2 ln(l + v^)) < exp(-Ci2t'/^S2). 

(4.37) 

The integral (13. 7p with fo{v) = fo,E/2{v) is finite by the same argument 
as for the moderate deviations. The extra factor (14.371) together with 
S2 > {si + S2)/2 leads to the bound exp(— (si + S2)) for t large enough. 

Case o"i > (72. Using the representation as in the moderate deviation 
case, we have, with respect to (Xi = e/2, the extra factor 

exp(itai ln(l - ^Ij2)) < exp(-ci3t'/'si), (4.38) 

from which we get the bound exp(— (si + S2)) as before. □ 

Proposition 7 (Uniform convergence on bounded sets). Fix any L > and 

Xi,Si with the above rescaling. Then, uniformly for (si,S2) G [—L,L]'^, 

lim(t/2)i/32^=^-^i /a;i-X2-l\ 
t^oo yn2 — n\ — \) 

= , ^ expf- ^/;~'^n l(r2>ri). (4.39) 

v/47r(r2-ri) V ^{r2-r^)) 

Proof. It is a special case of the first part of Proposition 5.1 of [1], where p 
is chosen such that k = 2~^/^ and (r^, Si) are replaced by (r,, Si). □ 

Proposition 8. For any Si, S2 G IR and T2 — ti > fixed, the bound 

(^/2)i/32-2-xi f^i -X2-l\ ^ ci2e'^''-''^ (4.40) 
\n2 -Tii-lJ 

holds for t large enough and c\2 independent oft. 

Proof. It is a special case of the first part of Proposition 5.5 of [1], where p 
is chosen such that n = 2~^/^ and (rj, s,) are replaced by (r^, s,). □ 
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5 About the other three regions 



(1) Constant density region. 

To obtain the result in tlie constant density region we consider the scaling 

X, = [{l-a)t/2-2n{t/2f/^-s,{t/2f'^] (5.1) 

with a > 1 fixed. The rescaled and conjugate kernel is as before 

irr^(ri, Si; r2, S2) = Kt{m,xi- n^, X2){t/2Y/'2"K (5.2) 

The binomial term is easily estimated and controlled. The main term Kl^^'^ 
is given by the formula (13.71) . with Si = Si, fi, f2, fs as in (13.81) . and the new 
/o is 

foiv) = -v + ln(l + v)-- \^{-v). (5.3) 

The two critical points of /o are now distinct, namely 

V- = —a/2 < —1/2 = v+. The constraint between the integration paths 
— 1 — Fq C r_i can not be satisfied if we want to choose Fq and r_i opti- 
mally, i.e., passing by and f_ respectively. For the analysis, one considers 
another representation of K^'^^'^, the same used in (I4.27p . The first term is as 
before but with the constraint F_i C — 1 — Fq and the second is the residue 
at M = — 1 — 1>, namely equal to I2 in (I4.28p . 

The first term is now controlled by choosing optimal paths for fo{v) and 
— /o(m), which pass by w+ and w_ respectively. fo{v+) < /o(f_), thus the 
first term is of order C(exp(t/o(u+) — t/o(f-))) = C(e^"*) for some a > 0. In 
particular, for a > 2 the first term vanishes identically (for t large enough), 
and as a \ 1, the first term is 0{e~^^°'~^^ /^^). 

The second term is just I2, up to some 2^^'^ factors due to the slightly 
different rescaling, the same kernel appearing in (5.5) of [2], where we already 
proved the pointwise convergence. The moderate and large deviations are the 
ones of I2 in (14.281) analyzed in Propositions and El In the t ^ 00 limit 
one then obtains 

hm Knn,si-r2,S2) = 2~'/'K^,i2-'/'n,2~'/'s,-2-'/'T2,2-'/'s2) (5.4) 

»oo 

with Kj\^^ is the kernel of the Airyi process. 
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(2) Linearly decreasing density region. 

To obtain the result in the linearly decreasing density region we consider the 
scaling 



m = [at/4 + r,(t/2)2/3], 



(5.5) 



with < q; < 1 fixed. The rescaled and conjugate kernel is 



(1 - v^/2)-2- 



■n2-xx-nx 



The main term of the kernel Kl^^'^ writes as fl3.7p with Sj = Sj — jo?!'^^ 
as in (13.81) . and 

fo{v) = -t; + (l-v^ + a/4)ln(l + t;)-(a/4)ln(-w), 
h{v) = (l-2/v^)ln(l + t;)-ln(-t;) + ln(v^/2) 

-(1-2/v^) ln(l- v^/2), 
f2{v) = -ln(l + t;) + ln(l- v^/2). (5.7) 

The function fo{v) has a double critical point at f = — \/a/2. The factor 
l + u + v = l — a/o at the critical point and the paths Fq and r_i can be 
chosen such that 1 + u + v remains uniformly bounded away from 0. The 
leading term of the integral comes from the neighborhood of the critical point. 
There, one applies the following change of variables, 

^ = -— + — ^iAW^ — " = + — ¥iwr^ — ^- (^-^^ 

Set Sh = a-2/3(2 - ^y^/^ and = a^/\2 - v^)-^/^. Then, the leading 
term in the main term of the kernel becomes 



(27ri)2 J J U -V e^V3+ri5feC/2-siS„c/ ' 
Thus 

KI''''''{ti, Si; T2, S2) SyK^^^Shn, S^si, ShT2, 8^82) (5.10) 
as t — >■ 00. By adequate control for moderate and large deviations, one proves 
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(3) Finite distance from the right-most particle. 

From the discussion on the initial condition, in particular from (12.121) . it 
follows that the asymptotic result is unchanged if one considers step initial 
conditions instead of our initial conditions. In [12] the case of step initial 
conditions was analyzed in a closely related model (a kind of discrete time 
TASEP but from the growth point of view). For step initial conditions, we 
have 

Kt{nx,Xx\n2,X2) = -\ l[„2>„,i + i^:4(ni, xi; n2, X2) (5.11) 

with 

n2 — 1 
fc=0 

where 

n{^) = j^^^fu{x^1n^), <l>r(y) = G(2/ + 2n,t), (5.13) 

the Cfc being the Charlier orthogonal polynomials [15]. This is obtained in 
the same way as in Appendix B of [2]. '^'^i^z) is the same as (B.7) of [2] with 
z — k replaced by 2; = x + 2n, and consequently the matrix S'fc,/ becomes the 
identity matrix. 

The Charlier polynomials converge to the Hermite polynomials H]. as 
follows 

lim(2t)'=/2^fc(t- V2ta,t) = {-\fEk{-cy)=}lk{a). (5.14) 
The scaling we have to use is 

Ui, Xj = [t - \p2ts^ (5.15) 

and the kernel rescaled as 

-52/2+52/2 

^r''(^l> Si; n2, S2) = V2t ^„,/2-ni/2 ^1' ^2, X2). (5.16) 



It is easy to see that the binomial contribution converges to 

S2-s,r--^-%s,>s,]. (5.17) 



g(s2-s2)/22(n2-ni)/2 



(ns - Til - 1)! 
Also, by fl5.14p . we have 



k/2 



hm *^(x.) = —-—=H,is.), hm $^(x.) = ^±-^H,{s.). (5.18) 
(2t)'=/2v27rt k\ 
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The kernel is a finite sum, thus 



t— »oo 

j=-n2 



where hk{s) = TT-^l^k\~^l^2~^/^Hk{s). 
(15T7D plus (15A9D gives 



(5.19) 



lim Kr\ni, si; ns, S2) = K'-'"'{n2, S2; ni, si) (5.20) 

with ii'^^UE ^YiQ kernel defined in Definition 1.2 of [12]. (Here we just order 
the entries differently). 

A Explicit form of the limit kernel 

Transition kernel in terms of Airy functions 

Let us denote 

Si = Si-mm{0,Tiy, Si = Si + max{0,TiY. (A.l) 

Then 

-^^00(7-1, si; T2, S2) = Ko{ti, Si; T2, S2) + Ki{ti, Si, T2, S2) + -^'2(^1, Si, T2, S2) 

(A.2) 

where 

Ko{n,Si;T2,S2) = — /. , exp -— — l[r,>r,h 

ggT^+ri.! ^47r(r2 - Ti) V 4(r2 - n) J 

(A.3) 

poo 

K2in,si;T2,S2)= dAe^("^-"i)Ai(s2 + A)Ai(Si + A), (A.4) 
Jo 

and 

poo 

Ki{ti,si;t2,S2)= / dAe^("2+"i)Ai(s2 + A)Ai(si-A). (A.5) 
Jo 

Equivalently, one can see that 

/■° 

Ki{n, si; T2, S2) = - dAe^("^+"^)Ai(s2 + A)Ai(si - A) 

^ — 00 

+2-1/3 (2-1/3(5^ + 5^ + _ ^2)2)) g-i(n+r2)(s2-5i) ^A.6) 
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B Trace class and transition kernel 



Proposition 9. The Fredholm determinant 

det(l-Xs^ooXs)w, (B.l) 

with Koo given in \2. 7| j or M.^i -z TO]) . is well defined, because there exists a 
conjugate kernel ofXsKooXs which is trace-class onTi = L^({ri, . . . , TmjxlR). 

Proof. In this proof, let us choose a Tq such that — Tq < ti < T2 < ■ ■ ■ < 
Tm < Tq. Denote by K'^°'^^ a conjugate of XsKooXs- Let Pk be the projector 
onto the space {/ G n\f{n,x) = for I =^ k} and Kl°^^ = PkK^°'^^Pi. Then 



m 



< 2^ II^M ii- (B.2 

From ( ]A.2I) we have 



ra=0,l,2 



where the conjugation function p(r, x) 7^ will be specified later. 

The formula defining KQ{Tk,x;Ti,y) is particularly nice in the 
Sk = Sk — [t^]^ variables (with [x]- = x for x < and otherwise). Thus 
we use the variables Sk instead of Sk- This is just a shift of the coordinate 
at the corresponding "time" r^. Thus, if we prove that Kl°P{Tk,x;Ti,y) = 

Kl°l^\Tk,x-\- [rfc]?.; Ti,y+ [t^]^) is trace-class for all Sk, k = 1, . . . ,m, bounded 

from below, then K^"^-' will also be trace-class for all Sk bounded from below. 

Therefore, we now work with the kernels and choose the conjugation 
functions (p(rfc,x) = p{Tk,x + [t^]?^)) to be 

p{Tk, x) = (1 + xy^e^'^'^y. (B.4) 

We analyze separately the three parts of the kernel. Let Kn{Tk,x;ri,y) = 
Kn{Tk,x+ [Tk]'i]n,y+ [ti]1), n = 0,1,2. 

Part a) KQ{Tk,x;Ti,y). We have 

lh>..l / , {y-xf\{l+x^r 
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2\2l 



In Lemma A. 2 of [1], we proved that the operator with above kernel is trace- 
class on L^(]R). (Recall that t; > if and only ii I > k). 

Part b) K2{Tk,x;Ti,y). We have 

^2(Tfc,a;;T,,|/)4^1[x>.-,]l[,>.-,]= / dXA,{x, X)A2{X,y) (B.6) 
p[n,y) Ju 

with 

Ai{x, A) = l[^>5,]l[A>o]P(Tfe, x)e-^>'^Ai{x + A + r^) (B.7) 
and ^ 

^2(A,?/) = l[A>o]l[^>3-,]^z^^Ai(y + A + rf). (B.8) 

Then we use IIA1A2II1 < ||^i||2 ||^2||2- Thus we have just to prove that 
Ai and A2 are Hilbert-Schmidt operators. This is easy to see, since 

||Ai||^ = / dxdX\Ai{x,X)\'^ (B.9) 

roo roo 

^ dx dXp{Tk,xfe-^^>'^\Ai{x + X + Tl 

Jsk Jo 

< C(To,5fc)<oo 



.2m2 



because the integrand is bounded, and for large x and A the decay is super- 
exponential due to the Airy function (Ai(2;) ~ e~3^ for z ^ 1). Similarly 
one shows that ||^2||2 < 00. 

Part c) Ki{Tk,x,Ti,y). We have 

i?i(Tfe,x,r,,y)^|;^l[,>,-,]l[,>,-,] = / dXB^{x,X)B2{\y) (B.IO) 

with 

B^{x, A) = l[.>.-,]l[A>o]ei"^'(l + x'')'"'e^>'^e'^°^Ai{x + X + r^) (B.ll) 

and 

B2iX,y) = l[,>ojl[^>,,je-^'^^-J-^/(A)^(A,|/) (B.12) 

with /(A) = e^n+Tk-2To)x ^nd ^(A, y) = e-^'ye-^°^Ai{y - A + rf). 

We need some estimates now. Since Ti + Tk — 2To < 2t^ — 2To, we have 

|/(A)|<e-^^ (B.13) 
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for fi = 2 (To — Tm) > 0. Moreover, 



Setting z = y + Ti 




(B.14) 



^(A,i/)|<cie--'^e-^«^|Ai(;^-A)|. 



(B.15) 



The first case is z < \. There, |Ai(2; — A)| < 1, thus 



(B.16) 



The second case is 2; > A (recall that A > 0). There 



^7(A,y)|<cie^«(^-^)Ai(z-A)<C3 



(B.17) 



because maXa;>o e"^''^Ai(a;) = C3 < cxd due to the super-exponential decay of 
Ai(x) for large x. Thus by (IB. 161) and (1B.17I) we conclude that, for all A > 
and y > si, there exists a constant C4 such that \g{X,y)\ < C4. 

The inequality II-B1II2 < 00 is similar to the ||v4i||2 case (use the decay 
of the Airy function). To see that II-B2II2 < 00, we use the bound flB.lSp to 
control the behavior in A, \g\ is just bounded by a constant and the decay in 
y is controlled by the (1 + y^)"^' term. 

In parts a), b) and c) we proved that all the kernel elements are trace-class 
on L^(]R) and this ends the proof of Proposition [HI □ 
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